From Luttinger liquid to Mott insulator: the correct low-energy description of the 
one-dimensional Hubbard model by an unbiased variational approach 
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We show that a particular class of variational wave functions reproduces the low-energy properties 
of the Hubbard model in one dimension. Our approach generalizes to finite on-site Coulomb repul- 
sion the fully-projected wave function proposed by Hellberg and Mele [Phys. Rev. Lett. 67, 2080 
(1991)] for describing the Luttinger-liquid behavior of the doped t—J model. Within our approach, 
the long-range Jastrow factor emerges from a careful minimization of the energy, without assuming 
any parametric form for the long-distance tail. Specifically, in the conducting phase of the Hubbard 
model at finite hole doping, we obtain the correct power-law behavior of the correlation functions, 
with the exponents predicted by the Tomonaga-Luttinger theory. By decreasing the doping, the 
insulating phase is reached with a continuous change of the small-q part of the Jastrow factor. 
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I. INTRODUCTION 

Electrons confined in one-dimensional (ID) systems ex- 
hibit peculiar non-Fermi liquid properties, that are by 
now rather well understoodi 2 ^ Let us just briefly men- 
tion some of them which will be useful for what fol- 
lows. Due to phase-space limitations, particle-hole ex- 
citations in single-band ID models are exhausted by col- 
lective charge and spin modes, which are dynamically 
independent, realizing what is commonly referred to as 
spin-charge separation. When these modes are gapless, 
they propagate as acoustic waves (zero-sounds), hence 
can be identified by two parameters, the sound velocity 
Ui and a dimensionless stiffness K^, i = p and i = a for 
charge and spin modes, respectively. Sometimes however 
the interaction opens a gap either in one of these sectors, 
for instance the charge sector in the Hubbard model at 
half-filling, or in both of them, like in a spontaneously 
dimerized chain. 

Besides spin-charge separation, another manifestation 
of the ID non-Fermi liquid behavior is the power-law de- 
cay with anomalous exponents of all correlation func- 
tions, when both spin and charge modes are gapless, 
or just some of them, when one of the two modes is 
gapped. Although these exponents are generically non 
universal, they all can be expressed in terms of the above 
mentioned K p and K a . Therefore the finite set of pa- 
rameters Up, Ucr, K p and K a is sufficient to characterize 
completely the asymptotic low-energy behavior of ID in- 
teracting electron models, similarly to the finite number 
of parameters which identify the low-energy behavior of 
Landau-Fermi liquids in higher dimensions^ Indeed, just 
in analogy with Fermi liquids, this kind of ID universal 
behavior was named "Luttinger liquid" by Haldanei£*§*£ 

In the case of non interacting electrons K p = K a = 1. 
If spin SU(2) symmetry is unbroken and the spin modes 
are gapless, K a remains one as for free fermions even in 
the presence of interaction and only K p parametrizes the 



anomalous exponents. In particular K p is smaller than 
one for repulsive interaction, and greater than one other- 
wise. When both charge and spin sectors are gapless, the 
asymptotic expressions of the charge and spin equal-time 
correlation functions are, apart from possible logarithmic 
corrections, 
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where n(x) and S(x) are the charge and spin density 
operators at position x, kp is the Fermi momentum, A\, 
Ai, and B are model-dependent constants. Analogously 
the singlet (and triplet) pairing correlations behave as 
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where A'(x) creates a singlet (or triplet) pair at position 
x. Finally, the non-Fermi liquid character of ID inter- 
acting models shows up transparently in the momentum 
distribution function near the Fermi momentum: 



-sign(fc — kp)\k — kp\ 



(4) 



where 9 is again expressed in terms of K p through the 
relation 9 = [K p — K^ 1 — 2)/4. For any finite interac- 
tion K p ^ 1, hence the momentum distribution function 
has a power-law singularity at the Fermi level, in con- 
trast to the finite jump characteristic of Fermi liquids. It 
is worth mentioning that, for non- frustrated models with 
repulsive interaction, the insulating phase at half-filling is 
characterized by all correlation functions decaying expo- 
nentially to zero apart from the spin-spin ones, which still 
decay as a power-law, formally like in J2J with K p = 0. 
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The actual value of the parameter K p depends upon the 
particular microscopic model, through the form and the 
strength of the interaction as well as the electron doping. 
For simple Hamiltonians K p can be explicitly calculated 
by the Bethe ansata^ or by exact diagonalization on small 
systems £ 

Although analytical techniques, especially bosoniza- 
tion£, give important insights into the low-energy proper- 
ties of ID systems, both in gapless and in gapped phases, 
yet they do not provide a simple representation of the 
ground-state wave function (WF). Even in those simple 
models which are solvable by Bethe ansatz, the actual 
ground state WF turns out to be very difficult to inter- 
pret and to deal with. In the Hubbard model, for in- 
stance, the ground state WF is very involved within the 
Bethe ansatz formalism and only in the strong-coupling 
limit it is possible to obtain significant simplifications be- 
cause of the explicit factorization of the WF into a charge 
and a spin part. 10 For this reason there have been many 
attempts aimed to find out approximate variational WFs, 
which, from one side, could correctly reproduce the pe- 
culiar properties of ID systems but, from the other side, 
could also give a transparent interpretation of their phys- 
ical properties>ii*i2ii2ii4iiii For instance it has been shown 
that the strong-coupling limit of the Hubbard model, 
where the on-site Coulomb repulsion U prohibits doubly 
occupied sites, is well described by a projected Slater de- 
terminant^ In particular the large-C/ conducting phase 
at finite hole doping has been found to be properly rep- 
resented through a fully-projected Gutzwillcr WF sup- 
plemented by a long-range density-density Jastrow fac- 
tor. The latter is a crucial ingredient which allows to 
recover the anomalous power-law behavior of the corre- 
lation functions. Indeed any short-range Jastrow factor 
cannot reproduce the correct Luttinger-liquid behavior 
since it is unable to affect the low-energy physics. For 
instance, it is well known that the simple Gutzwiller WF 
has a finite jump in the momentum distribution func- 
tion^ which is certainly incorrect in ID. 

In this paper we generalize this approach to a finite 
Coulomb repulsion U, namely to the case where charge 
fluctuations are still allowed. We show that even at fi- 
nite U it is possible to design a consistent WF, which 
can faithfully describe the evolution from the Luttinger- 
liquid behavior at finite hole doping to the Mott insu- 
lating phase at half-filling. Again the crucial ingredient 
turns out to be a density-density Jastrow factor applied 
to a simple Slater determinant. More specifically, we ap- 
ply this variational WF to analyze the single-band Hub- 
bard model. Thanks to the recent improvements in the 
energy minimization schemes for correlated WFsjii the 
long-range tail of the Jastrow factor can be determined 
very accurately, without imposing any parametric form, 
even if the change of its tail may contribute to a very tiny 
energy gain. We will show that the low-energy properties 
of the Hubbard model are correctly reproduced within 
this unbiased variational approach. 

The paper is organized as follow: in Sec. [H] we intro- 
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FIG. 1: The value of (3 from Eq. JTTJ for half-filling and 
quarter-filling for different ratios U/t. 



duce the model and the physical quantities we are inter- 
ested in and in Sec. lIIII we present and discuss the results. 



II. THE MODEL 

As we previously mentioned, in this work we aim to 
design a variational WF which is equally accurate for 
both weak and strong Coulomb repulsion. For sake of 
simplicity, we will test the quality of this WF in the well- 
known single-band Hubbard model on a finite chain with 
L sites, N particles and periodic boundary conditions. 
The Hamiltonian is 
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where we use the standard notations in which Cj l<T (c\ a ) 
destroys (creates) an electron with spin a at the site i 

and m,a = 4,a c hv 

The variational ansatz for the ground state WF is 



\^n) = J\V), 



(6) 



where J is a Jastrow factor and I'D) an uncorrelated 
Slater determinant, that, for simplicity, we assume to 
be the iV-electron Fermi sea of the tight-binding model 
with dispersion efe = — 2icosfc. Since spin SU(2) is un- 
broken, there is no need to introduce a spin-spin Jastrow 
factor, because a free-electron determinant already pro- 
vides the correct value K a = 1. Therefore, in order to 
preserve all the symmetries of the Hamiltonian, J is a 
purely density-density correlator written in term of the 
density operator m = J2a c l a c i-° as 
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FIG. 2: Charge structure factor for L — 82 and U/t = 10 FIG. 3: Left panel: charge structure factor N q for the Hub- 

at half-filling: comparison between the variational result and bard model at quarter-filling for L — 124 and different U/t. 

the RPA expression given by Eq. Hill . Inset: zoom around Right panel: the same for the spin structure factor S q . 
q = 7T. 



where, due to translation and inversion symmetry, v% t $ 
can be described by using L/2 independent variational 
parameters v(\i — j\) = v^y, in particular, v(L/2) can 
be set to zero since, due to the conservation of the total 
number of particles, Vij — > Vij + const provides only an 
irrelevant normalization factor in the WF. 

A more sophisticated possibility would have been to 
choose \V) as the ground state of a mean-field BCS 
Hamiltonian projected onto a fixed number of particles. 
For the simple Hubbard model of Eq. J5J), the latter 
ansatz leads to a slightly lower energy, without modifying 
the low-energy properties, and, therefore, it will not be 
considered in the following. It should be mentioned, how- 
ever, that in presence of a large next-nearest-neighbor 
hopping term, a gap in the BCS mean field Hamiltonian 
may open, leading to translational symmetry breaking at 
half-filling or to a phase of singlet pairs with dominant 
superconducting correlations at finite doping. 18-19 

Summarizing, all the variational parameters are con- 
tained in the Jastrow coefficients v(\i — j\) for all the L/2 
independent distances in real space and are calculated by 
a full energy minimization, without assuming any partic- 
ular parametric form. 

Important quantities to assess the nature of the ground 
state are the static structure factor for the charge: 
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and similarly for the spin: 
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where n q and S q are the Fourier transform of the local 
density and spin operators. 

Another quantity that gives information on the WF is 
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FIG. 4: Log-log plot of the quasiparticle weight Zk at k = 
kF — vr/4 as a function of L in the quarter-filled Hubbard 
model for different values of U/t. Lines are power-law fits. 



the quasiparticle weight: 



Zk = 



Kttjv-i|c fc , g |gjy)| 2 

;*iv|*jv)(*iv-i|*jv-i>' 



(9) 



where |*jv) and |*jv-i) are the WFs with N and (JV— 1) 
particles, Ck,a is the annihilation operator of a particle of 
momentum k and spin a. The (JV — l)-particle state is 
obtained by the JV-particle one by removing an electron 
from the Slater determinant, i.e., I^at-i) = J Ck,a 
In a Fermi liquid Z k is finite in the thermodynamic limit 
signaling the existence of coherent quasiparticles. On the 
contrary a non-Fermi liquid phase without quasiparticles 
is identified by a vanishing Z k . 
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III. RESULTS AND DISCUSSION 



The important role of the Jastrow factor has been al- 
ready discussed by Hellberg and Mele in the context of 
the ID t—J modelii^ In this case, it is possible to show 
analytically that the momentum distribution function of 
the variational WF has an algebraic singularity at k F , 
with an exponent related to the strength of the Jastrow 
factor iSi Recently, we have showni^ that in the ID Hub- 
bard model a suitable long-range Jastrow can also drive 
a metallic Slater determinant into an insulating phase 
and that, in the limit of small momenta, the Fourier 
transform of the Jastrow parameters v q are related to 
the charge-structure factor. In particular, in the con- 
ducting regime and and for small momenta a random- 
phase-approximation (RPA) expression holdsi^i 



N° 
g 



1 + 2v q N° ' 



(10) 



where we indicate with N® the charge structure factor 
of the WF without J . For a conducting Slater determi- 
nant, like the free-electron state, N q = \q\/Tr for small q. 
For a non-vanishing hole doping S and interaction U/t, 
by optimizing the variational WF, we obtain that the 
Jastrow factor is singular, i.e., v q ~ This fact is 

crucial to recover the correct low-g behavior of N q , whose 
linear slope is renormalized by the interaction, leading to 
N q ~ Kp\q\/n. Therefore, the Jastrow factor has to be 
intrinsically long-range. According to the RPA expres- 
sion UlUfl . a more singular Jastrow term v q ~ 1/q 2 is 
needed to bring the system into the insulating phase at 
half-filling^ In ID, such a Jastrow factor leads to a con- 
fined phase, where perturbation theory does not apply, 
and, therefore, also the RPA is expected to fail. Nonethe- 
less, it turns out that, also within the insulating phase, 
the small-<7 behavior of N(q) is qualitatively reproduced 
by Eq. I|1U|I and a more accurate empirical expression is 
given by: 
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where (3 > 2 is a constant that strongly depends upon the 
electronic interaction. In Fig. ^ we report (3 for different 
values of the ratio U/t at half-filling, and, for comparison, 
also at quarter-filling, where the value (3 — 2 is recovered, 
according to Eq. (jTU|) . In the half-filled insulating phase, 
although for g^Owe have that N q ~ q 2 , the coefficient 
of the quadratic term is not simply related to the Jastrow 
factor. In this case, it is important to emphasize that the 
presence of the singular Jastrow factor v q ~ l/<7 2 deter- 
mines a qualitative change of the static structure factor 
N(q) even at large momenta q ~ 2kp = 7r. In fact, the 
charge structure factor N°(q) for the free Fermi gas is 
characterized by a cusp at q ~ 2kp, which is responsi- 
ble of the well-known Friedel oscillations in a metal and 
the RPA expression cannot remove this feature. On the 
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FIG. 5: Log-log plot of the pairing correlation function at the 
maximum distance P(L/2) as a function of L in the quarter- 
filled Hubbard model for different values of U/t. Lines are 
power-law fits. 



other hand, a true insulating phase does not possess this 
singularity and, indeed, as shown in Fig. [21 the charge 
structure factor N(q) for the WF containing a singular 
Jastrow factor shows a smooth behavior around q ~ ir. 
This clearly indicates the non-perturbative and highly 
non-trivial effect implied by the formation of a confined 
state in this correlated WF between empty sites (holons) 
and doubly occupied ones (doublons). 

In order to demonstrate that the WF © is able to 
capture the Luttinger-liquid metallic properties, we con- 
sider the quarter-filling case. In Fig. |3| we show the 
charge and spin structure factor for different values of 
U /t. For small momenta, the linear slope of N q is renor- 
malized with respect to the non-interacting value, lead- 
ing to N q ~ K p \q\/-K. On the other hand, the small-g 
behavior of S q is not affected by the interaction and we 
have that S q ~ |g|/47r. Notice that, in the presence of a 
strong interaction, the two singularities at 2kp and 
are clearly visible in N q , whereas in S q , only the singu- 
larity at 2kF can be detected. From the small-g linear 
part of N q , it is possible to extract the value of K p (see 
Table P), which is in very good agreement with the exact 
one£ It is important to emphasize that the relationships 
among exponents of different correlation functions are 
correctly reproduced by our variational WF. Indeed, we 
can compare the value of the exponent 6 found from a di- 
rect evaluation of the quasiparticle weig ht © at k = k F , 
i.e., Zk ~ 1/L e (see Fig. @}, with the one obtained with 
9 = (K p + K~ 1 — 2)/4 by using the value of K p extracted 
from the linear slope of N q . As reported in Table []] wc 
obtain an excellent agreement for the values of the inter- 
action U/t considered. Finally, we can also calculate the 
singlet pairing correlations 



P(r) = i£(#iv|A i+r A 



(12) 
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TABLE I: Critical exponents for the ID Hubbard model at 
quarter- filling: K p is found from the low-g behavior of N q , 
9c = (K p + Kp 1 - 2)/4, and 9 is found by fitting Z k with 
Zk ~ 1/L e . The last two columns refer to the critical expo- 
nent of the pairing correlations: a is found from the pairing 
correlation at the maximum distance P(L/2) ~ 1/L a and 
a c — K^ 1 + 1. In the first column, we report the exact value 
aiK p . 

U/t K e p xa " K p 9 9 C a a c = 

~ 4 0.711 0.705(3) 0.031(5) 0.031(3) 2.1(1) 2.42(6) 
10 0.594 0.595(3) 0.078(5) 0.072(3) 2.4(1) 2.68(9) 
18 0.551 0.550(3) 0.097(5) 0.092(3) 2.5(2) 2.82(9) 
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FIG. 6: Left panel: charge structure factor N q for L = 82, 
U/t — 10, and different hole dopings. Right panel: N q for 
5 — 2/11 and different L. 

where 



creates a singlet pair of electrons at nearest neighbors. 
In order to calculate the exponent a related to the de- 
cay of P(r) ~ l/r a , we consider the pairing correlation 
at the maximum distance P(L/2) for different sizes, see 
Fig. [SJ In this case, the signal is very small and a precise 
determination of the critical exponent is quite difficult. 
Nonetheless, the results reported in Table [I] are rather 
satisfactory and not too far from the ones obtained with 
the exact relation a — K~ + lA 

Let us now consider how the insulating case is reached 
by decreasing the hole concentration. First of all, it 
should be mentioned that, not too close to the insulating 
phase at half-filling, the charge and spin structure factor 
have small size effects and, therefore, reliable calculations 
are possible even without using too large L. As an exam- 
ple, we report in Fig. (right panel), the case of doping 
5 = 2/ 1 1 , where we can see that there are no appreciable 
differences in N q from L = 22 to L = 1 54. 

In the doped region, the system is always conducting, 
N q having a linear behavior for small momenta, with 
a slope that depends upon U/t and 5. For sufficiently 



FIG. 7: Jastrow factor v q (multiplied by q) for different hole 
dopings, obtained by a careful minimization of the energy for 
U/t = 10 and L = 82. 

small hole doping, it turns out that the linear regime 
is limited to a small window around q — 0, whereas for 
larger momenta, N q acquires a finite curvature, see Fig. El 
(left panel). The two different regimes are separated by 
the singularity at q = Akp = 2nS, and, therefore, by 
decreasing 5, the width of the linear regime shrinks, the 
slope being almost constant. Therefore, we arrive at the 
empirical result: 

N q ~ ^le(4fc F - q) + (c+ q 2 )Q{q - 4fc F ), (14) 

where kp = (1 — 6)w/2, 0(x) indicates the Heaviside 
step function and c is a constant determined by imposing 
continuity of N q at q = 4fcp. This singular behavior, with 
the kink at q — 4fc F , is entirely due to correlation and it is 
compatible with the exact result that K p remains finite, 
more precisely K p — > 1/2, for 5 — -> 0.- It is important to 
stress that, by approaching half-filling, the width of the 
linear regime shrinks with doping and reduces to zero for 
6^0. 

For completeness, we report in Fig. the form of the 
Jastrow factor at half-filling and for small concentration 
of holes considered in Fig. El Starting from the insulating 
phase, upon doping, v q moves away from v q ~ 1/q 2 , and 
becomes less singular, i.e., v q ~ l/\q\. Clearly, at very 
small doping, the size effects affect the small-g part, and, 
in particular, for the smallest momentum we can have 
some deviations from the expected v q ~ l/|<z| behavior. 

The correct minimization of the Jastrow factor is par- 
ticularly important for having an accurate description of 
Zk, especially when approaching half- filling. Indeed, in 
this case, the Jastrow factor for one hole is considerably 
different from the insulating one for small q's (see Fig. |HJ) , 
and one has to optimize both WFs with N and (N — 1 ) 
particles^ This difference can be appreciated by con- 
sidering = l/L J2 q ^o w ?' w hich diverges linearly with 
the system size L if v q ~ 1/q 2 and, instead, diverges only 
logarithmically if v q ~ 1/q. In Fig. |8J we report Si as 
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FIG. 8: Jastrow factor v q for one hole, U/t = 18 and 
different L. Lines connect the points Inset: the quantity 
S_t = l/^E,^o v i ( see text) as a function of L for half-filling 
and one hole. Lines are fits of the data. 
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FIG. 9: Log- log plot of Zk as a function of L for U/t = 10 
and 18. Lines are power-law fits. 

a function of L for the insulating state and for the one- 
hole case: the difference between the two cases clearly 
demonstrates the different behavior of v q for small mo- 
menta. By a careful minimization of both the WFs, it is 
possible to recover the result that 6=1/2 independently 



of U/t. 2 ?- Indeed, upon increasing U/t, our variational 
WF gives a rather accurate description of the insulating 
phase, the size effects being strongly reduced due to the 
small correlation length expected at large U/t. In this 
limit, we obtain a reasonable good agreement with the 
exact exponent for the quasiparticle weight (see Fig. [§J| : 
9 = 0.60 ± 0.05 and 9 = 0.55 ± 0.05 for U/t = 10 and 
U/t = 18, respectively. On the other hand, it should be 
mentioned that a naive calculation with a singular Jas- 
trow v q ~ 1/q 2 for both WFs would lead to a wrong 
exponential behavior of the quasiparticle weight. 

In conclusion, we have shown that an accurate descrip- 
tion of the charge and spin properties of the ID Hubbard 
model is possible within the simple variational ansatz 
given by Eq. JJjJ. The main ingredient of the variational 
WF is represented by a singular density-density Jastrow 
factor, whose long-range part determines the low-energy 
behavior of the correlation functions and can be found 
by optimization of the energy. In the conducting regime, 
the small-q Fourier transform of the density-density Jas- 
trow term diverges like v q ~ implying the correct 
renormalization of the charge structure factor for small 
momenta, i.e., N q ~ K p \q\/iv. By approaching the insu- 
lating phase and by decreasing the hole doping, the Jas- 
trow factor modifies its small-g part and, eventually, it 
acquires a more singular behavior like v q ~ l/q 2 ■ Within 
this context, it is possible to obtain a consistent sce- 
nario, where we recover the correct relationships among 
the exponents of different physical quantities. Obviously, 
our approach works due to the quality of the variational 
ansatz in ID. Despite the difficulty to define very accu- 
rate WFs in higher spatial dimensions, the possibility to 
obtain the correct low-energy description of a Hamilto- 
nian by using a simple optimization of the energy is very 
appealing and work is in progress for the generalization of 
this ansatz to two-dimensional systems, where no exact 
solutions are available. 
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